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The fast decay of carrier inversion in photoexcited graphene has been attributed to optical phonon 
emission and Auger recombination. Plasmon emission provides another pathway that, as we show 
here, drives the carrier relaxation dynamics on ultrafast timescales. In studying the nonequilibrium 
relaxation dynamics we find that plasmon emission effectively converts inversion into hot carriers, 
whose energy is then extracted by optical phonon emission. This mechanism not only explains the 
observed fs-lifetime of inversion but also offers the prospect for atomically thin ultrafast plasmon 
emitters. 
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Graphene owes its extraordinary electronic, optical 
and plasmonic properties [2-6] to the presence of Dirac 
points in its band structure, in the vicinity of which elec¬ 
trons are described as massless Dirac fermions (MDFs). 
One of graphene’s unique characteristics is that low- 
energy pair excitations can decay into plasmons [7, 8], 
a process that is kinematically forbidden in most 3D and 
2D materials. 

When excited with a femtosecond optical pulse, a hot 
carrier plasma is generated, which rapidly thermalizes 
into a state of quasi-equilibrium due to carrier-carrier 
scattering on 10-fs scales [9, 10]. After thermalization, 
electrons and holes recombine via interband scattering 
processes, effectively depleting carrier inversion over time 
[11, 12]. Femtosecond pump-probe experiments [IS¬ 
IS] and time- and angle-resolved photo-emission spec¬ 
troscopy (tr-ARPES) [16-18] indicate the presence of two 
dominant recombination channels: a slow decay on a 1-ps 
scale due to emission of optical phonons [19-21], whose 
signatures can be observed by Raman spectroscopy, and 
a fast decay on a 100-fs scale that has been attributed to 
Auger recombination (AR) [15, 22]. This interpretation 
remains subject to discussion as AR processes are su- 
pressed when considering dynamic screening in random- 
phase approximation (RPA) [23, 24]. 

In this Letter we show that nonequilibrium plasmon 
emission (NPE) (see Eig. 1) provides an alternative path¬ 
way for carrier relaxation that drives the decay of inver¬ 
sion on a 100-fs scale. It was recently established that 
plasmon emission in inverted graphene is ultrafast, with 
rates that exceed those of optical phonon emission by 
at least one order of magnitude [25, 26]. However, in 
assuming predefined carrier and plasmon distributions, 
these works did not account for the dynamic evolution of 
the coupled carrier-plasmon system far from equilibrium. 
Thus, despite the prediction of ultrafast rates, it remains 
an open question how NPE impacts the relaxation dy¬ 



FIG. 1. (Color online) Schematic of carrier recombination 
channels in inverted graphene, (a) Nonequilibrium plasmon 
emission and (b) Auger recombination both convert elec¬ 
tron/hole pairs into hot carriers. In contrast to Auger re¬ 
combination, plasmon emission is an inelastic process that is 
enabled by the strong coupling of electron/hole pair to 2D 
plasmon excitations. 


namics, as plasmons are not only constantly emitted but 
also absorbed back into the electron/hole plasma. Ulti¬ 
mately, to assess the importance of NPE as channel for 
carrier relaxation, it is necessary to study the temporal 
evolution of the nonequilibrium plasmon distribution in 
interaction with the inverted carrier system, and together 
with optical phonon emission, which provides additional 
pathways for carrier recombination and cooling. As we 
will show, an interplay of channels emerges, where NPE 
drives the ultrafast conversion of inversion into hot car¬ 
riers on 100-fs scales, while optical phonons efficiently 
extract heat from the MDE plasma. 

As basis for this work we introduce nonequilibrium 
rate equations for an inverted MDE plasma in interac¬ 
tion with bosonic reservoirs via emission and absorption 
processes. In quasi-equilibrium, the state of the inverted 
carrier plasma is characterized by chemical potentials 
{a = e, h) for the electrons and holes and a com¬ 
mon plasma temperature Tc, whose temporal changes 
can be expressed by those of the carrier densities = 
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N{fi^,Tc) of electron and holes and the total energy den¬ 
sity U = ^ oi the plasma [where = U (/i^, Tc)] 

(see [ 1 ], Sec. I). We split the time-derivatives into optical 
excitation (pump) and relaxation terms, 

= ^Ipump - ^Irel , and ^7 = f7|pump - Ulrel , (1) 

where as particles and holes are created (and 

annihilated) pairwise. For quasi-monochromatic excita¬ 
tion with frequency huj and intensity Iuj{t) the pump 
terms are given by A/'lpump = Ke[crinter(^)]|^u;(t)|V( 2 ^^) 
and L^lpump = Ro[cr(co’)] |F^c^(t)p/ 2 . Here Zq is the vac¬ 
uum impedance, = 2ZQ\t{uj)\‘^I^^{t) the electric 

field square, t{uj) = (1 -f Zoa{uj)/2)~^ the transmission 
coefficient, a = a{uj; ^Tc) the optical conductivity, 

and (Jinter its interband only contribution. 

Carrier relaxation, in turn, is a result of the interaction 
with bosonic reservoirs, such as the plasmon and optical 
phonon fields [27]. The microscopic mechanisms behind 
the dynamic relaxation of the carrier number and energy 
densities are intra- and interband scattering processes, 
as given by Boltzmann collision integrals. To account for 
the various bosonic reservoirs (z/), as well as intra- and 
intraband transitions (A), we write 

NU = J2R^,eh, U\rel = Yl E ’ ( 2 ) 

v' V X=ee,hh,eh 

The net carrier and energy relaxation rates Rj^^eh and 
are obtained by summation over all wavevector 
states q of the reservoir bosons, i.e., 

^ E ^ E > (3) 

q q 

where 

»'r^,A(q) = 7tA(q)[^>^(q) + M- (4) 

is the net spectral emission rate, A the area of the 
graphene sample, the energy-dispersion 

of the boson field, 7 ^^^ ilux) emission (absorption) 
rate, and ^^^(q) the nonequilibrium distribution of the re¬ 
spective boson field. The latter evolves in time according 
to 

n^(q)= E '’>^.A(q) - nr.(q)ldecay (5) 

X=ee,hh,eh 

with a decay term that accounts for additional loss chan¬ 
nels. Within the relaxation time approximation we de- 
fine n,.(q)ldecay = intro- 

ducing the decay rate and the equilibrium distribu¬ 
tion n^f^^(e,To) of the reservoir at ambient temperature 
Tq. Note that, as ^^^(q) returns to equilibrium, r^^xi^) 
becomes zero according to the condition of detailed bal¬ 
ance. 



FIG. 2. (Color online) (a) Plasmon frequency dispersion 
ujp\{q) and (b) loss spectrum 7 pi(ct;) of intrinsic graphene 
for varying temperature Tc scaled with chemical potential fi. 
The dispersion crosses through regimes of interband gain (red 
area) and absorption (blue area). For /x = 0.1 eV tempera¬ 
tures shown range from Tc = 0 — 3208 K. 

Equipped with these equations, we proceed to study 
the carrier relaxation dynamics for interaction with the 
plasmon reservoir and then under the inclusion of phonon 
channels. The results in this work are obtained by in¬ 
tegrating Eqs. ( 2 ) and (5) over time. Eor all our sim¬ 
ulations we consider a suspended sheet of graphene at 
ambient temperature that is excitated by a pulse with 
a fiuence of 133/iJ/cm^ and a photon energy of leV. 
To enforce the same initial conditions for all simulations, 
we keep all relaxation channels switched off during op¬ 
tical excitation. Under these conditions the excitation 
pulse creates a thermalized inverted state with chemical 
potentials ~ 0.3 eV and a carrier temperature 

Tc 2288 K. 

When the carrier system of graphene is inverted, elec¬ 
tron/hole pair excitations and plasmons interact strongly 
via plasmon emission and absorption processes. The as¬ 
sociated spectral rates can be calculated approximatively 
using Eermi golden rule (EGR), with an accuracy that 
critically depends on the exactness of the plasmon fre¬ 
quency dispersion ujp\{q) [26]. Therefore we trace the 
exact complex-frequency roots uj{q) = ujp\{q) — ijp\{q) of 
the dynamic dielectric function, i.e., the solutions of 

eRPA(g,w) = 1 - = 0. (6) 

Here Vq = e^/(26og) is the bare 2 D Coulomb potential 
and ^^{q,uj) the irreducible nonequilibrium polariz¬ 
ability o^ inverted graphene at temperature Tc. Eigure 2 
shows frequency dispersion ujp\{q) and the loss spectrum 
7 pi(<^) = 7 pi(^pi(<^)) foi* temperatures kBTcl{2fi) in the 
range of 0 — 1.5. Plasmons experience interband gain for 
hwpi < 2 /x (red area) and inter- and intraband loss for 
hwpi > 2 /x (blue and grey areas) at all temperatures, as 
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FIG. 3. (Color online) Relaxation dynamics of the coupled 
electron/hole and plasmon system, (a) Chemical potential 
fi (black) and carrier temperature Tc (red) for TcoII = 100 fs 
(solid lines) and varying collision time (shaded area), (b) 
Temporal evolution of nonequilibrium plasmon distributions 
np\{q) and emission/absorption edge (black dashed lines). 


evident from the change of sign of the net stimulated ab¬ 
sorption rate 7pV™(<^) = 27 pi(co’). Note, for our analysis, 
it is sufficient to consider interband processes only, as 
reabsorption above the Fermi-edge rapidly depletes the 
plasmon mode population before entering the intraband 
regime. The respective interband emission/absorption 
rates are given by 


7pl,eft(9) 


2ag0{uj - vpq) 



^Re[£RPA(g,u;)] 

duj 


■ (7) 

uj=ujpi{q) 


Here ag = afc/vp ~ 300/137 is the fine-structure con¬ 
stant of graphene and (g, cj) a dimensionless 

measure for the phase-space of the absorption and emis¬ 
sion processes (see [1], Sec. II). In addition to these pro¬ 
cesses, one needs to consider the loss of plasmons due to 
collisions [28], which emerges as a result of various elas¬ 
tic and inelastic velocity scattering processes [29] . Values 
for Tcoii are typically in the 50-500 fs range, depending 
on carrier temperature, chemical potentials and impu¬ 
rity concentration. While collision loss only weakly im¬ 
pacts on the frequency dispersion cjpi(g) [26], it serves as 
a secondary decay channel as plasmons purged from the 
reservoir are not available for reabsorption at later times. 

We first study the interaction of carriers with plasmons 
in isolation, i.e., without inclusion of the optical phonon 
channels. Figure 3 shows the evolution of the chemi¬ 


cal potential (/r = /i® = /i^) and the carrier temperature 
Tc for varying collision times, together with the temporal 
traces of the nonequilibrium plasmon distribution npi(g). 
Within the first 100 fs after photoexcitation, a burst of 
plasmon emission occurs that leads to a drop of inver¬ 
sion (solid black line) to roughly half its initial value; at 
the same time, the carrier temperature (solid red line) 
increases by around 500 K [see Fig. 3(a)] as plasmons are 
emitted below the Fermi-edge. Over the next 200 fs the 
recombination process gradually slows down. This is be¬ 
cause reabsorption of plasmons above the Fermi-edge and 
the decrease of the plasmon emission rate with temper¬ 
ature creates a plasmon emission bottleneck. The tem¬ 
poral traces of the nonequilibrium plasmon distribution 
npi(g,t) [Fig. 3(b)] are shown together with the thresh¬ 
old wavevector gth(^) (black dashed line) for which the 
net plasmon gain is zero, i.e., 7pi(gth) = 0 (no plasmon 
gain/loss). Below this threshold plasmons are emitted, 
above they are reabsorbed. Collision loss partially re¬ 
moves the emission bottleneck as it purges plasmons from 
the reservoir and thus prevents reabsorption of plasmons 
into the electron/hole plasma. As a result, the recom¬ 
bination of carriers accelerates with increasing collision 
rate [grey shaded area in Fig. 3(a)]. 

Having analyzed the plasmon channel in isolation, we 
next study the interplay of NPE and optical phonon emis¬ 
sion, which plays a pivotal role as they facilitate the re¬ 
combination of electron/hole pairs [20], the cooling of 
hot carriers [21], as well as collision loss [30]. In in¬ 
verted graphene, or at high temperature, the inter- and 
intraband emission of longitudinal and transverse optical 
(LO/TO) phonons are dominant channels for carrier re¬ 
combination and cooling. For this work we consider all 
relevant optical phonon channels, the FO, the KO, and 
the KA phonons with energies of epo = 196 meV, cko = 
160meV, and cka = 120meV, respectively [31, 32]. 
As the optical phonon modes are quasi dispersion-free, 
Eq. (3) reduces to Rjy^x = + 1] — and 

5'zy,A = where is the occupation number and 

Cy the phonon energy. Closed-form expressions for the 
rates F^^ = and the phonon density of states 

My^x are given in the supplement [1], Sec. HI. 

To understand the fundamental difference of NPE and 
optical phonon emission, we first consider the carrier 
dynamics without NPE, i.e., under inclusion of optical 
phonon emission only assuming a phonon decay time of 
Tiat = 2.5ps [21]. Eigure 4(a) depicts the temporal evolu¬ 
tion of the chemical potential (black line) together with 
the carrier temperature (red line) and temperatures of 
the FO, KO phonons (blue lines; KA phonons omitted 
for clarity). The initial dip in carrier temperature (at 
t ^ 400fs), accompanied by a rise of phonon temper¬ 
atures and chemical potential, is predominantly due to 
intraband phonon emission, which continuously extracts 
energy from the carrier plasma. As the phonon tempera¬ 
ture equilibrates with the plasma temperature, intraband 
























4 




FIG. 4. (Color online) Relaxation dynamics of carriers coupled to (a) phonon only; (b) plasmons and phonons where TcoII = 
100 fs. Shown are chemical potential // (black line), carrier temperature Tc (red line), and temperatures of the FO, KO phonons 
Tro (blue solid line), and Tko (blue dashed line); the grey area in (b) indicates change with collision time. (c) + (d) Effective 
decay times Ti and amplitudes Ai extracted from tri-exponential fit to /i(t) considering plasmons only (dashed lines), plasmons 
and phonons (solid lines); grey lines indicate decay times extracted from (a) (phonons only). 


cooling becomes increasingly inefficient. Over the next 1 - 
ps, /i drops from 0.3 eV to 0.2 eV mainly due to interband 
emission of optical phonons. As carrier temperature and 
chemical potential further decrease both inter- and in¬ 
traband emission of optical phonons slows down. Less 
electron/hole pairs are available at the required phonon 
energies and carrier cooling is bottlenecked by the slow 
decay of optical phonons into acoustic phonons on ps- 
scales. 

The relaxation dynamics changes dramatically when 
NPE is taken into account [see Fig. 4(b)]. The rapid 
drop of inversion within the first 100 -fs, due to the plas- 
mon emission burst, is followed by a gradual slowdown as 
plasmons are reabsorbed above the Fermi-edge. The plas- 
mon energy that flows back into the electron/hole sys¬ 
tem heats the carrier plasma and thus prevents the drop 
in carrier temperature that was observed in Fig. 4(a). 
At t = 1 ps the chemical potential has fallen well below 
0.1 eV, compared to a value of /i ^ 0.25eV in Fig. 4(a), 
where plasmon emission was switched off. Collision loss 
further accelerates the decay of inversion, albeit not as 
strongly as in Fig. 3(a) (plasmons only), as phonons now 
provide an efficient cooling channel that alleviates the im¬ 
pact of carrier heating due to plasmon reabsorption. The 
combination of plasmon and optical phonon emission ef¬ 
fectively bypasses bottlenecks observed for the isolated 
channels, thereby accelerating the carrier recombination 
and cooling dynamics. 

To analyze the decay of carrier inversion due to 
NPE we fit /r(t) with tri-exponential functions /ifit(t) = 
Mo The extracted (effective) decay 

times Ti and relative amplitudes Ai are shown in Fig. 4(c, 
d) in dependence on collision time. NPE alone induces 
three timescales (dashed lines): a fast rate (ri ~ 30fs) 
that relates to the plasmon emission into unoccupied 


plasmon modes with an amplitude that rises quickly for 
^coii < Ti; a slower rate (ri ^ 300 fs) due to fast emis¬ 
sion and reabsorption at the Fermi-edge; and a slow rate 
(ts ~ 3ps) that is strongly influenced by absorption be¬ 
low the Fermi-edge, as evident from the rise of ampli¬ 
tude As with Tcoii- Activating the phonon channel in 
addition to NPE (solid lines) causes an interplay of chan¬ 
nels as the NPE decay times r 2 and ts are close to the 
timescales that govern the decay via phonon emission 
[Fig. 4(c); grey lines]. This is apparent in the change 
of amplitudes A 2 and A 3 , which for TcoIi > 200 fs be¬ 
come almost constant, resulting in a decay of inversion 
that is almost independent of TcoIi [see Fig. 4(b); grey 
shaded area]. Most prominently, the relaxation dynam¬ 
ics and extracted timescales show that inversion above 
Me = = 0.1 eV decays on 100 -fs scales due to the 

ultrafast decay of electron/hole pairs into plasmons. 

In conclusion, we have established that plasmon emis¬ 
sion drives the ultrafast carrier relaxation in photo- 
excited graphene, as 2D plasmons can couple strongly 
to pair excitations of the inverted carrier plasma. Our 
results are consistent with the recent experimental obser¬ 
vation of fs-decay of population inversion [14, 16]. In con¬ 
trast to Auger processes, whose experimental detection 
is very challenging [18], plasmon emission can be directly 
observed and thus provides a novel path for experimen¬ 
tal characterization of relaxation processes. Interaction 
of plasmons with a (polar) substrate causes a red-shift of 
the plasmon emission spectrum and hybridization with 
surface optical phonons. These effects are not expected 
to change the main conclusions drawn here, and will be 
studied in a future work. 
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